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Inflationary reheating is almost completely controlled by the Floquet indices, Hk- Using spectral 
theory we demonstrate that the stability bands (where /i^ — 0) of the Mathieu and Lame equations 
are destroyed even in Minkowski spacetime, leaving a fractal Cantor set or a measure zero set of 
stable modes in the cases where the inflaton evolves in an almost-periodic or stochastic manner 
respectively. These two types of potential model the expected multi-field and quantum backreaction 
effects during reheating. 
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I. INTRODUCTION 

Inflation is perhaps the heir apparent in high energy 
cosmology. No other theory can be cast successfully in so 
many guises. Yet, to be successful, any inflationary sce- 
nario must reheat the universe. This epoch is typically 
very short and has become a powerful virtual laboratory 
for researching non-equilibrium quantum field theory in 
curved spacetime [p]-pT| . It is defined by rapidly evolving 
paradigms (see the left half of Fig. 1), attempting to deal 
with the non-perturbative processes occurring in preheat- 
ing . The defining characteristic of a preheating model 
are its Floquet indices ^fc, which control the number of 
produced particles rik ^ [ dkk^e 
the fields |9|. 

Of particular interest are the issues of how the /i^ vary 
as (i) the parameters (typically the coupling constants) 
of the theory are changed, and (ii) the functional form 
of the infiaton evolution is altered. The first has been 
studied in depth in the case where the infiaton evolution 
is exactly periodic and is the very origin of preheating. 
The second issue is much less studied and is the subject of 
this paper where we introduce spectral methods as an el- 
egant classification tool for inflationary reheating classes; 
though it also applies to the parametric amplification of 
gravitational waves To illustrate the idea, consider 
the Hill equation: 



J rffcfc^e^™'^''* and variances of 



ij+[A- 2qP{2t)]y = 



(1) 



where P{2t) is any periodic function of the independent 
variable t. When P(-) cx cos(-), we have the Mathieu 
equation, while P(-) = cn^(-) yields the Lame equation, 
where cn(-) is the elliptic cosine function. Now Eq. (Q) is 
equivalent to the one-dimensional Schrodinger operator- 
eigenvalue problem: 



dx 



-JZ2 + Qi^)y = 



(2) 



under the transformations x t, Q{x) 2qP{2t), A <-!■ 
A. We denote the spectrum of C by (j(C), its complement 
by cr(£). The crucial point to notice is that the set of 
modes with positive Floquet index of Eq. (|^) contains 



the set (7{C) of Eq. (g). By exploiting this equivalence 
we will show that: 

I - Generically for almost-periodic potentials Q{x), the 
corresponding spectrum is a fractal Cantor set. This im- 
plies a Floquet index which is positive on a dense subset 
of momenta k S [0, oo). As the strength of the couphng 
is varied, drastic changes to the spectrum can occur that 
are absent in the purely periodic case. 

II - Preheating exists in stochastic infiation and the 
Floquet index is positive with Lcbesgue measure 1 for all 
momenta k. Estimates for fik are given in the perturba- 
tive and non-perturbative limits by Eq.'s (^5|) and ( p^ ) 
respectively. 

These results describe geometrically the breakup of the 
stability-instability chart that in the case of the Hill equa- 
tion (|^) has a neat band-structure. Visually, this is sim- 
ilar to the breakup of invariant tori in the KAM theory 
of chaos . Similar results to those of II above have re- 
cently been found independently in specific cases 
for the noise. In particular, Zanchin et al [ pT| rewrite Eq. 
(|l|) as a first order matrix equation and then use Fursten- 
berg's theorem regarding products of independent, 
identically distributed random matrices to obtain esti- 
mates of , showing that noise increases the /ifc over the 
periodic case and that fj,k > for all wavelengths, thus 
also demonstrating the break-up of the stability bands, 
although in the context of small noise. 

The results of I and II above enlarge the known re- 
heating territory (see Fig. 1). The classical theory of 
reheating was developed in the early 80 's with recog- 
nition of the importance of perturbative resonances p3[ . 
In this case the resonance bands essentially dwindle into 
a discrete sequence of lines. The paradigm shift of pre- 
heating 01, showed that the resonance bands could be 
very broad at large q, a fact that extends to gravita- 
tional waves was followed by approxima- 
tions to deal with the non-perturbative effects ||,^ , and 
the idea of non-thermal symmetry restoration induced by 
the large quantum fluctuations Recently stochastic 
resonance HJl^ has represented a shift away from the 
simple picture of static resonance bands, due to the large 
phase fluctuations which, mod 2tt, behave randomly in 
the non-perturbative regime in an expanding universe. 
Finally, the negative-coupling instability B represents 
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the exciting possibility that modes can evolve where the 
spectrum of the corresponding Schrodinger operator is al- 
most empty ^ , and Floquet indices are very large. This 
mechanism has been realised within the framework of 
non-minimal fields coupled only to gravity, which there- 
fore offers a geometric reheating channel due to the 
oscillation of the Ricci curvature. 
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FIG. 1. A schematic map of models and approximations in 
preheating. Minima[ references for techniques in the exact[y 
periodic case are shown in braclcets. The right-hand branch 
corresponds to paradigms deve[oped in this paper. 

In this paper we consider almost-periodic (section ^) 
and random (section |l|) potentials which arc both real- 
izations of metrically transitive (i.e. ergodic and homo- 
geneous) operators [po|,pl|. 



II. CANTOR REHEATING 

We limit our discussion in this section to Minkowski 
spacetime |^. Consider the inflaton, 0, and the two min- 
imally coupled scalar fields ip, x with the natural mass 
hierarchy ^ rrix- In this section we consider 

the effective potential: 

Vi^,^,x)-^<t>' + -^^' + ^^' + ^x' (3) 

F{-),J{-) are assumed to be analytic in their respective 
arguments. The evolution of the Xk modes is then M: 



Xk+\^-^+m'^+g'Fi(b) + h'Jiip)jXk = (4) 

Instead the inflaton zero-mode evolves according to: 

4> + ml^f^^ + X(b'' ^0 (5) 

where the frequency of oscillation is partially controlled 
by the effective mass: 



*This can be restated as the fact that the integrated density 
of states, proportional to the rotation number [M, vanishes. 

^The general expanding FLRW case is contained within the 
spectral theory of Sturm-Liouville operators. 



'.eff 



F' 



(6) 



In the case that (x^) = {5<t>^) = 0, the inflaton simply 
oscillates with constant period. However, when there are 
multiple- fields, or the mass acquires corrections due to 
quantum fluctuations, the period is no-longer constant, 
but may increase monotonically, oscillate or exhibit ran- 
dom fluctuations. As an example, in the special case of a 
Yukawa interaction, F{(f>) = 0, Eq. (H) has a pure driving 
term cx {x^)- 

Understanding the behaviour of the Floquet indices in 
these more general cases is extremely complex and we 
therefore search for special, solvable classes, which nat- 
urally leads to the study of almost-periodic potentials 
p5[ . In this case, the spectral theory becomes signifi- 
cantly richer and more beautiful. 

In general, for a separable Hilbert space, as we have 
here, the spectrum can be decomposed w.r.t. an abstract 
measure d/i as a = uac U crp U ctsc @]i where uac is 
the absolutely continuous, ap is the pure point and asc 
the singular continuous part of the spectrum. An impor- 
tant constraint is that the Floquet index must vanish for 
A G (JAC, he. GAc = {A e R I /iA = 0} bar a set of 
measure zero [Q, and hence aAc corresponds to modes 
with bounded evolution (extended eigenfunctions). On 
the complement, ctac, the Floquet indices are known to 
be positive, and hence in the case when oac is empty, 
fj,k > almost everywhere p9[ |. 

To understand this better, consider the periodic poten- 
tial. In this case, a{C) = uac and on <j{C), (X\ > 0. It is 
the fact that aAC is so large that forces the strong band 
structure o f the Mathieu and Lame equations. In the 
Lame case, cr(£) consists of a single band and there are 
very few exponentially growing modes jj] . A comparison 
between the archetypal Cantor set and the spectrum of 
the Mathieu equation is shown in Fig. (||). 

The almost-periodic potentials arc those for which 
the Fourier transform consists of a frequency basis {wi} 
where the smallest vector space containing this basis, Ai, 
is dense in R. If M is generated by finitely many uji, 
then the potential, Q{t), is quasi-periodic, i.e. Q{t) = 
f{ujit, ...,cj„t) with /(<i +mi, ...,t„ -|-TO„) = f{ti, --^n) 
with TTii £ Z and the oji pairwise incommensurate [ p5[ . 

As a simple example consider Eq. (^) with F{(l)) = 
(jy^ , J{(p) = and A = 0. Then the equations for the 
quantum fluctuations of Xk are: 



Xk ■ 



fc2 
«2 



9 



2 J,2 



(7) 



with (f> ~ sm{m^t) and ^ siii{m^t). When m^/rriip is 
irrational, the potential is quasi-periodic. The spectrum 
of Eq. (g) for almost-periodic potentials can be pure 
point and is, in a non-rigorous way, generally a nowhere 
dense Cantor set Hence, in the example above, for 
infinitely many irrational values of m^/m^, we may ex- 
pect that the spectrum of Xk will be nowhere dense, and 
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consequently that almost all modes will grow exponen- 
tially. 

However, unlike the random potentials to be consid- 
ered in section { f]% , the Lebesgue measure of the spec- 
trum, even if a Cantor set, need not be zero |25| and aAc 
may be non-empty. Indeed it is possible to have a — uac 

H. 
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FIG. 2. A schematic diagram comparing the spectrum 
of the Mathieu equation (upper) and that of a generic 
limit-periodic operator with fractal Cantor spectrum (lower). 
The black shading represents the complement of the spec- 
trum, i.e. modes with /i^ 7^ 0. The Cantor spectrum is only 
shown at second level due to resolution limitations. The same 
pattern repeats itself in each white gap and hence the number 
of stable modes is much smaller than in the periodic case. 

More rigorous results exist in the case that Q is limit- 
periodic, i.e. it is a uniform limit of periodic potentials. 
A typical example is provided by: 



00 

Q{t) a„ cos 

n=l 
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< 00 



(8) 



For these potentials we can call on rigorous theorems: 

Theorem 1 - a{C) of Eq. is generically ^ a 

nowhere dense Cantor set for Q an element of the space 
of limit periodic potentials. Hence ct{C) is dense in R. 

Thus the set of k for which Hk 7^ 0, is dense in R+ 
despite the fact that aAc is not empty. This existence 
of only a Cantor set of stable modes leads us to call this 
cantor reheating. An associated issue is what happens 
as the coupling to the potential (analogous to q in Eq. 
(1)) is increased. While in Eq. (1) this simply changes 
the breadth of the instability bands and magnitude of 
fik while always leaving a = o^ac, in the almost-periodic 
case this is not the case. Indeed, as in the example be- 
low, the nature of the spectrum (w.r.t. the splitting into 
cp, o'sci <^Ac) can change suddenly with q. 



ing studies [|6[|j]. A special case is the almost-Mathieu 
equation with exhibits a rich variety of effects ||: 

— y{x -|- 1) — y(x — 1) + 2qcos{ax + uj)y{x) — Xy{x) 

(9) 

When a is rational the potential is periodic and one has 
the band spectrum of the Mathieu equation. However, 
for irrational a, the spectrum is a Cantor set generically 
for pairs {q,a) £ [Q. This shows that the nowhere 
dense nature of the spectrum is not lost as one increases 
q and hence moves from perturbative reheating to broad- 
resonance preheating. Further, when a is irrational, the 
Floquet index has the lower bound fSSt]: 



IJ-kiq) > ln|g| 



(10) 



so that for \q\ > 1, fj,k > and hence the absolutely con- 
tinuous part of the spectrum, ctac, becomes empty. In 
this case, if a is a Liouville number [**|, then the spectrum 
is purely singular continuous, a = use ||2^- Conversely 
if Igl < 1, the point part of the spectrum is absent for 
irrational a |34l . 



III. STOCHASTIC INFLATIONARY REHEATING 

We now consider the case where the potential Q in 
Eq. (|^) is random. This models for example, the clas- 
sical limit of stochastic inflation, where the dynamics of 
the local order parameter in a FLRW background, are 



described by 



F2 



^ + iH<j, + V'{ci>)^—V"'{4>)m. 



(11) 



where H = a/a is the stochasticaly evolving Hubble con- 
stant ||2^ , ^ is a coloured gaussian noise of unit amplitude 
with a correlation time of order H~^, so that the infla- 
ton evolves stochastically. The origin of the noise is the 
backreaction of quantum fluctuations with wavelengths 
shorter than the coarse-graining scale ^3 24 p^. 

Here we will consider the potential V{(j)) — \(f)'^ /A and 
as before an interaction term g^cjp'x^ l^.. The quantum 
fluctuations of the fields (j) and x then given by 0: 



A. The discrete almost-Mathieu equation 

Consider the discretized Schrodinger eq. (H), which is 
also used (in higher dimensions) in numerical preheat- 



■'■The space of limit periodic potentials is a complete metric 
space and hence generic means here "for a dense G5" e.g. feg. 



^For example, Eq. (|9|) also exhibits the beautiful property of 
duality, similar to the S-duality of string theory, since under 
Fourier transform g — > 1/g and A — > —X/q 

**A Liouville number, a, is irrational but well approximated 
by rationals so that there exist integers Pn,qn — > 00 and a 
number C with \a — p„/qn\ < Cn~''" . 

^^Here we neglect the backreaction of the produced field x 
on the expansion; valid during the first phase of reheating 
before backreaction terminates the resonance. 
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fc2 
n2 



4-^-3 



(a 



,3/2; 



/'fe)=0, 

(12) 



^+m^ + ^p + 5VM(«^/V)-0. (13) 



where p — K{^<j)'^ — V{(f))) is the pressure and k — 8ttG. 
These equations are again equivalent to Eq. (^ , but this 
time with stochastic potentials, which are the opposite 
extreme to periodic potentials since they generically have 
empty ctac- In ^ct we have the following theorem: 

Theorem 2 - If Q{x) is a sufficiently random 

potential for Eq. (0), then fi\ > for almost all X E H 
and a AC of C is empty with probability one. 

Indeed we see that a positive Floquet exponent is guar- 
anteed for all modes bar a set of measure zero, and again 
reheating is significantly different from the periodic case. 
Here "sufficiently random" means nondeterministic ]2C| ] 
and is typically a requirement that correlations decay suf- 
ficiently rapidly. An example is a Gaussian random field 
whose correlation function has compact support. To pro- 
ceed to obtain quantitative estimates of the Floquet in- 
dices, we exploit the fact that Eq.s (p^ , [l3| ) have the form 
of stochastic harmonic oscillators. 



A. Explicit estimates for the Floquet indices 

Consider the stochastic harmonic oscillator, with fre- 
quency given by uj'^ = +q^{t), where g is a dimension- 
less coupling to the mean-zero coloured stochastic pro- 
cess ^{t) and — jc? in our case. The Floquet index 
has been shown to be strictly positive for all q It 
has also been related to the spectral density of fluctua- 
tions via averaging over the second moments ||3^ of the 
random process ^ pq]: 



Cmat-t'))cos[2{uj^)t']dt' 



(14) 



This has the form of a fluctuation-dissipation theorem 
since fluctuations in the inflaton field determine the dis- 
sipation rate into other fields. In the case that £,{t) is a 
mean zero ergodic Markov process, fik can be explicitly 



estimated as |36 



0{q' 



(15) 



where / is the Fourier transform of the expectation value 
of the two-time correlation function {^{t)^{t — t')). Note 
that ^jLk oc so that although all modes grow expo- 
nentially, the Floquet index is, as in the periodic and 
stochastic resonance |^ cases, a rapidly decreasing func- 
tion of fc. In the broad-resonance limit, g ^ oo, we write 
= Kq + qnl, which gives (assuming max ^ < k^) js^: 



47r 



de / d^- 



0(i/Vg) 



G(ln(K2 -^cos^e)) 
(16) 



where G is the infinitesimal generator of ^(t) defined by 
the limit of the operator sequence: G — limt^o{Ut — 
T)/t. Here I is the identity operator and Ut is a family of 
operators on the space of bounded continuous functions 
/, defined by Utf{x) = E[f{mMO) = where E 
denotes the expectation operator |2^. Again /i^ > for 
all fc, and given a model of evolution, one can explicitly 
estimate /ifc , and hence the numbers of produced particles 
Hk and variances {S(f>'^) and (x^)- 



IV. AMPLIFICATION OF GRAVITATIONAL 
WAVES 



Here it is demonstrated that the evolution equatons 
for gravitons can be cast in the form of the Schrodinger 
Eq. (|^). Within the Bardeen formalism ||l^], or by us- 
ing the Weyl tensor, there exists a strong correspondence 
between gravitational waves and scalar fields during re- 
heating |14|. In the Bardeen approach, gravitational 
waves are described by transverse-traceless metric per- 
turbations, with mode functions hk satisfying [ pT| : 

ia'^'h,r+(^^ + lpyay^hk)^0, (17) 

so that Q — 3p/4, k'^/a^ ^ A and a'^^^hk y es- 
tablishes the equivalence with Eq. (^) . This implies that 
during reheating with an almost-periodic or stochastic in- 
flaton evolution, all wavelengths will be amplified, rather 
than just those in the Floquet instability bands. In the 
stochastic case with small noise, Eq. ( |l5|) shows that the 
tensor spectrum will acquire a tilting with p,k ~ fc~^, 
causing further subtleties for the inflationary potential 
reconstruction program. 



V. CONCLUDING REMARKS 

Using spectral theory we have shown the breakup of 
the stability bands of the Mathieu equation. For constant 
q, the set of stable modes becomes a Cantor set generi- 
cally for almost-periodic inflaton evolution and a set of 
measure zero if the evolution is sufRciently random. This 
implies that quantum corrections to the effective mass 
will lead to an exponential growth of all modes of the 
reheated field, and of the gravitational wave background 
produced during inflation. Further, there exists the pos- 
sibility, absent in the periodic case, of drastic qualita- 
tive changes to the nature of the spectrum at certain 
transitional values of q, as demonstrated by the discrete 
almost-Mathieu equation. 

This is of relevance to the issue of non-thermal sym- 
metry restoration, since the quantity Xt(£-) f-kdk is much 
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larger in the classes considered here. This should re- 
flect directly in the number of particles produced, rifc, 
the non-thermality of the spectrum and the variances of 
the fields. 

We have largely limited ourselves here to the study of 
spectral changes as the potential is changed and found 
drastic changes from the simple band structure of Flo- 
quet theory. However, as a warning, beautiful counter- 
examples exist, though they are suitably rare. Firstly, 
using inverse scattering theory for Eq. (g), one can show 
that there exist, and infact derive all examples of, almost- 
periodic potentials which have band spectra similar to 
the periodic case. An example is given by the Bargmann 
potentials ||3^ . Hence it is not always true that the spec- 
trum of an almost-periodic operator must be a Cantor 
set. 

Neither is it true that altering the coupling to the po- 
tential always induces a change in the spectrum. From 
the theory of isospectral deformations of the Schrodinger 
equation it is known that if Q{x,y) is a potential of Eq. 
(g) and satisfies the KdV equation Qy—dQQx+Qxxx = 0, 
then (t(£) is independent of y and there thus exist 1- 
parameter families of potentials with the same spectra 
pO| . The connections with integrable systems become 
clear, and hence so do the rarity of such events. 
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